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Pade-approximant methods are used to extract information about leading positive zeros or poles of QCD 
and SQCD j3- functions from the known terms of their perturbative series. For QCD, such methods are seen to 
corroborate the flavour-threshold behaviour obtained via lattice approaches for the occurrence of infrared-stable 
fixed points. All possible Pade-approximant versions of the known (one- to four-loop) terms of the QCD MS 
)3- function series are consistent with this threshold occurring at or above n/ = 6. This conclusion continues to be 
true even if higher-degree Pade-approximants are introduced to accommodate an arbitrary five-loop contribution 
to the QCD /3-function for a given number of flavours. 

The dynamics characterising the infrared region of asymptotically-free gauge theories are becoming an area 
of active investigation. Historically speaking, the property of asymptotic freedom has been linked with "infrared 
slavery" at large distances - the non-observability of gauge-group nonsinglet particles as asymptotic free-particle 
states because of the growth of the gauge coupling constant to arbitrarily large values at some infrared momentum 
scale. Such methodological simplification, however, rests upon the existence of a Landau singularity in the evolution 
of the gauge couplant, a singularity that is likely an artefact of truncating the gauge couplant's /3-function to any 
given order. Consider, for example, the QCD /3-function to two-loop order for three active flavours: 
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Given some empirical initial value, e.g. x T — a s (m T )/7r , eq.(l) has the analytic solution 
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As [i decreases from m T in this expression, the couplant x(/i) grows from x T to become infinite at 

fi L = m T (1 + 9/(16x r )) 32/81 exp [-2/(9s T )] 
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the Landau singularity. The temptation is to identify /i^ with the infrared boundary of QCD as a perturbative 
gauge theory of quarks and gluons. Note that if a s (m T ) = ttx t — 0.35, the central PDG value [1], then [ij, is found 
from (3) to be 493 MeV. 

However, this Landau singularity is really a consequence of truncating the perturbative /3-function series (1) to 
its renormalization scheme independent two-loop order terms. Subsequent /3-function terms are renormalization- 
scheme-dependent and presumably negotiable (I will have more to say on this later on). For example, we can 
conjecture two toy-model /3-functions 

9 dx 9 x 2 
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which both agree to two- loop order with eq. (2), but have infrared behaviour manifestly different from the Landau 
pole characterising (2). The former equation has a /3-function pole at x = 9/16 = 0.5625. Consequently, as [i 
decreases from m T , x increases from x T to a maximum value of 0.5625, which is achieved at a critical infrared 
momentum scale 
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The domain of x(fi) is restricted to /i > fi c , in which case \x c represents an infrared-boundary of QCD at which 
the gauge couplant achieves its largest possible but finite value. The existence of an infrared boundary is no longer 
coupled with the arbitrary growth of the interaction coupling. 

By contrast, the /3-function of eq. (5) has a zero at Xirfp = — 8/(9i?2) + \/64/(81i?2) I/-R2 (R2 is negative), 
meaning that increases to eventually level off at xirfp as /i — > 0. For this case, which corresponds to the 
existence of an infrared-stable fixed point (IRFP), there is no infrared boundary to gauge-theoretical QCD - the 
domain of cc(/x) includes all (positive) values of fi. 

A point that needs to be made here is that if (1), (4), and (5) were to represent all-orders /3-functions arising 
from legitimate but differing renormalization schemes, [J then the infrared behaviour of the couplant would itself 
be scheme-dependent. Indeed, one could easily construct a /3-function candidate consistent with (1) to two-loop 
order whose higher order terms are sufficiently large and negative [e.g. an arbitrarily large negative R2 in (5)] 
to ensure that an IRFP occurs at an arbitrarily small value for x, suggesting that QCD in such a scheme would 
remain a tractable perturbative gauge theory at arbitrarily large distances. Since the infrared behaviour of QCD 
is empirically known to be confining at sufficiently large distances, such apparent scheme-dependence must be 
overstated. We therefore take the point of view here that differing but self- consistently realized renormalization 
schemes will necessarily lead to the same physical results - in particular, to equivalent infrared dynamics. 

A vector-like gauge theory that illustrates all this is that of N = 1 supersymmetric SU(3) Yang-Mills theory, 
supersymmctric gluodynamics characterised only by gluon and gluino fields. Because of its supersymmetry, the 
all-orders /3-function for this theory can be obtained algebraically by requiring that the anomaly multiplet within 
the theory respect the Adler-Bardeen theorem [3,4], or alternatively via instanton-calculus considerations [5]. Sur- 
prisingly, the same /3-function arises from either approach, and even more surprisingly this /3-function, like (4), is 
characterised by a pole: 

"«"M/i<* P) 

As noted in [6], such a pole necessarily implies a critical momentum scale /x c for the asymptotically free phase of 
the theory that would constitute an infrared boundary to the region for which the theory is perturbative. The pole 
value of the couplant [x(/j c ) = 1/6] constitutes an infrared-attractive point terminating the evolution of the (real) 
gauge coupling constant in its asymptotically-free phase. We henceforth denote (7) as the NSVZ /3-function, after 
the authors of ref. 5. 

The /3-function for supersymmetric gluodynamics has also been calculated to three non-leading orders within 
the dimensional-reduction (DRED) renormalization scheme [7]: 

= -9x 2 [l + 6x + 63x 2 + 918x 3 + ...] (8) 

The higher order terms of the /3-function (8) differ from those of the geometric series implicit in (7) once one 
gets beyond their equivalent leading and next-to-leading terms. Both NSVZ and DRED schemes are presumably 
valid ones (although only the former upholds the Adler-Bardeen Theorem to all orders). One can even obtain a 
perturbative road-map between couplants in the two schemes: if z = x NSVZ and y = x DRED , then [4] 

y = z [1 + az + (27 + 6cr + a 2 )z 2 + (351 + IVTa + 15a 2 + a 3 )z 3 + ...] (9) 

where the constant a is arbitrary as a consequence of both (asymptotically-free) couplants having identical leading 
and next-to-leading terms within their respective /3-function series]^] Consequently, unless one specifies initial values 
for couplants in each scheme, one cannot say anything at all about the relative size of couplants in the two different 
schemes at a given momentum scale. If, for example, the all-order extension of the DRED series (8) exhibits a pole, 
such a pole need not occur at its x = 1/6 NSVZ location. Nevertheless, consistent infrared dynamics does require 
that a pole indeed occur in the DRED scheme. In other words, the all-orders extension of (8) should not exhibit the 
Landau-pole or IRFP behaviour seen to characterise of the /3-function examples (1) and (5), respectively. Rather, 
the DRED scheme should the same infrared dynamics as are evident from (7): a finite infrared-attractive point 
x(fi c ) terminating the evolution of the asymptotically-free phase of the couplant, with [i c serving as an infrared 
bound on the domain of x(fi). 

1 Equation (1) the /3-function truncated after scheme-independent two-loop order terms, has already been used to generate an 
apparently self-consistent renormalization scheme [2]. 

2 An arbitrary constant in the relation between perturbative couplants in two different renormalization schemes of a given theory 
will occur provided the leading and next-to-leading terms of the /3-function series in both schemes are the same. Such is the case, of 
course, for conventional QCD. 
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To get some insight into the infrared behaviour of DRED supersymmetric gluodynamics, the four known terms 
of the /3-function series (8) have been utilised to construct Pade approximants - ratios of polynomials in x whose 
power-series expansions reproduce these known terms. Pade approximants are often employed not only to predict 
next order terms of a series, such as the (3- function series for QCD [8], scalar field theory [8,9,10], and supersymmetric 
QCD [4,11], but also to explore whether infinite series with only a few known terms can be expected to have zeros 
or poles [12]. |^] In reference [4], for example, it is demonstrated that every Pade approximant (except the truncated 
series itself) that reproduces the first four known terms of a series which differs infinitcsimally from the geometric 
series 1/(1 — |r|x) exhibits a pole that 

1) differs infinitesimally from the true pole at x — l/|r|, and 

2) occurs prior to any positive spurious Pade-approximant zeros. 

In other words, the pole driven infrared dynamics of the NSVZ /3-function (7) could have been predicted from Pade 
approximants constructed from that series' first four terms. 

Relevant Pade approximants whose power series reproduce the four known terms of the DRED /3-function (8) 
are [4] 

9 dx „ o (1 — 14x) 
u — n — -9x 2 - — i 10) 



2 dx n_2 (! ~ 8.5714x - 24.4286x 2 ) 
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dx „ « 1 
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The first positive pole of (10) occurs at x = 0.0604 and precedes the zero at x = 1/14. Similarly, the pole of (11) 
at x = 1/14.5714 = 0.0686 precedes that approximant's first positive zero (x = 0.0924), and the approximant (12) 
exhibits a first positive pole at x = 0.0773. Although (12) is constructed to have no zeros other than x = 0, there 
is no a priori reason that this approximant should exhibit a positive pole of comparable magnitude to those of (10) 
and (11). These results (and further analysis of higher approximants presented in [4]) are clearly indicative of the 
same pole-driven infrared dynamics that characterise the NSVZ /3-function (7) for the same N = 1 supersymmetric 
SU(3) Yang-Mills field theory. 

The techniques illustrated above are applicable to QCD itself. There is considerable controversy as to the 
infrared dynamics which characterise QCD, as well as the flavour dependence of such dynamics. The idea that 
the QCD couplant freezes out to some effective a s (fi = 0) at sufficiently low momentum scales appears to have 
both some theoretical justification [14] (even for the rif = case [15]) and as well as phenomenological utility [16]. 
However, such IRFP dynamics for rif ^ 6 are inconsistent with both theoretical studies based upon a /3-function 
series truncated after its first two scheme independent terms [17] as well as with a lattice study indicative of an 
rif = 7 threshold for IRFP dynamics [18]. 

A Pade-approximant approach to the infrared dynamics of QCD similar to the example presented above is 
formulated in detail in refs. [9], [13], and [19]. In this work, the known (and first unknown) terms of the QCD MS 
/3-function series [20] 

P 2 -^ = -f3 x 2 S(x), pa = 11/12 -n//6, x = a s {n)/^ (13) 



S(x) = 1 + [(51/8 - 19n//24)//3 ] x 

+ [(2857/2 - 5033n//18 -I- 325nJ/54)/64/3 ] x 2 

+ [(114.23 - 27.134/1/ + 1.5824^ + 5.8567 • 10" 3 n^)//3 ] x 3 

+ i? 4 x 4 + ... (14) 

are used to construct various Pade-approximants whose leading positive zeros/poles are compared to see which 
occur first (or occur at all) . We choose to use the MS renormalization shceme simply because the /3-function series 
in this scheme is known to higher order than in any other perturbative scheme - moreover, phenomenological QCD 
is overwhelmingly based upon MS calculations. We denote by S^ N ' M ^(x) the Pade-approximant to S(x) whose 

3 Straightforward examples of this method for recovering the first positive pole or zero of sec(x) ±tan(x) are presented in ref. [13]. 
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numerator and denominator are respectively degree-N and degree-M polynomials of the couplant x. If N + M = 3, 
the power-series expansion of the approximant S'[ JV I M ] (x) is of sufficiently high order to reproduce the known series 
terms in (14); e.g. for rif = 3 

cf2|il M _ l-2.9169x- 3.8750a; 2 

* lj ~ l-4.6947x ( ' 

S m (x) - 1 - 8.1734a: (16) 

b [) ~ 1-9.951 lx + 13.220^ Ubj 

In both of these approximants, a positive pole precedes any positive zeros, indicative of the same pole-driven 
infrared dynamics known to characterise NSVZ supersymmctric gluodynamics. If one constructs such approximants 
for any choice of rif (as in ref. [13]), one finds for rif < 5 that S 1 ^' 1 ' always has a positive pole preceding any 
positive zeros; moreover, the same statement applies to S^ 1 ' 2 ! as well, provided rif < 6. These results suggest that 
pole-driven infrared dynamics characterise QCD for up to five or six flavours. The occurrence of positive zeros that 
are not preceded by positive poles, (i.e., /3-function zeros corresponding to IRFP dynamics) does not characterise 
gPI 1 ] until rif > 7, and does not characterise S^ 1 ' 2 ' until rif > 9, corroborating the existence of a flavour threshold 
of IRFP behaviour [17,18] as well as the absence of such behaviour below this threshold [21]. 

In ref. [13], these same qualitative conclusions are shown to be upheld by S'[ JV I M 1 for N+M =5, i.e., for Padc- 
approximants to S(x) whose power series are constructed to reproduce the four known terms of (14) and the 
arbitrary five loop term R4X 4 . In such approximants, each numerator and denominator polynomial coefficient of 
powers of x is itself linear in the parameter R4 - e.g. for rif — 3 

[2 | 2] 1 + (7.1946 - 0.10261^4)2; + (-11.329 + 0.075644i? 4 )a; 2 

^ ~~ 1 + (5.4168 - 0.10261i? 4 )a; + (-25.430 + 0.25806.R 4 )a; 2 ( ' 

Regardless of the value R4 takes, one finds that positive zeros which precede every positive pole do not occur in 
the approximants S^ 3 \ S^ 2 ^ , or S^ 3 ' 1 ! until n/ > 9, 7, and 6, respectively, consistent with the non-occurrence 
of IRFP dynamics below these threshold values of rif. Moreover, for arbitrary R4 a positive pole which precedes 
every positive zero does occur in the approximants S^ 1 ' 3 ! and S^ 2 \ provided rif < 5. The approximant S^ 3 ' 1 ' has 
a single positive pole only if R4 > 0, a consequence of having a denominator linear in x, but this pole is found to 
precede any numerator zeros in the approximant for all positive values of R4 as long as rif < 7. 

These results are clearly indicative of the occurrence of pole-driven infrared dynamics for QCD once heavy 
flavours are decoupled. Similar pole-driven infrared dynamics are also shown in [13] to characterise QCD in the 
N c — > 00 't Hooft limit. Moreover, very recent work [18] has demonstrated that different approximants to the 
N c = 3, rif — 3 MS QCD /3-function exhibit surprising consistency in their predictions of infrared boundary 
coordinates (// c ,x(// c )) associated with pole-driven infrared dynamics. The implications of such dynamics, in 
particular the possibility of having both a strong and an asymptotically-free phase of QCD with common infrared 
properties [6], have only begun to be explored. 
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